Nanomaterial Characterization in Complex Media—Guidance and Application

A broad range of inorganic nanoparticles (NPs) and their dissolved ions possess a possible toxicological risk for human health and the environment. Reliable and robust measurements of dissolution effects may be influenced by the sample matrix, which challenges the analytical method of choice. In this study, CuO NPs were investigated in several dissolution experiments. Two analytical techniques (dynamic light scattering (DLS) and inductively-coupled plasma mass spectrometry (ICP-MS)) were used to characterize NPs (size distribution curves) time-dependently in different complex matrices (e.g., artificial lung lining fluids and cell culture media). The advantages and challenges of each analytical approach are evaluated and discussed. Additionally, a direct-injection single particle (DI sp)ICP-MS technique for assessing the size distribution curve of the dissolved particles was developed and evaluated. The DI technique provides a sensitive response even at low concentrations without any dilution of the complex sample matrix. These experiments were further enhanced with an automated data evaluation procedure to objectively distinguish between ionic and NP events. With this approach, a fast and reproducible determination of inorganic NPs and ionic backgrounds can be achieved. This study can serve as guidance when choosing the optimal analytical method for NP characterization and for the determination of the origin of an adverse effect in NP toxicity.

This development leads to safety concerns about this emerging group of functional materials [15][16][17]. The research in this area reveals that a broad range of parameters of the NPs have an influence on the outcome of NP-related toxicological studies. On the one hand, different intrinsic physico-chemical parameters of the NPs, such as size [18] and shape [19], as well as surface properties, such as charge [20] or chemical modifications [21], have an influence on the fate of the NPs and any successive biological response [22]. On the other hand, additional interactions of the NPs after their suspension, e.g., with specific media or artificial liquid constituents, may modify the NP behavior in the dispersion In this study, metal oxide NPs made of CuO were used as a model system due to their broad size distribution and pH-driven dissolution properties [42]. In combination with the model matrices, these samples were used to evaluate the capabilities of different techniques (DLS, microwave digestion with ICP-MS and direct-injection (DI) spICP-MS) in challenging scenarios. Measurements were conducted every hour for the first 8 h and then after 1, 2 and 7 days. The high frequency of time points at the start of the experiment takes into account the potential higher reactivity and aims to reduce possible errors based on these. A DI spICP-MS method was developed and optimized to directly measure small volumes of the sample without the requirement of additional preparation steps. The development process was based on the technical specifications ISO TS 19590:2017 [43]. The approach developed in this study can be used as a fast method to determine inorganic NPs and their potential health risks attributed to the dissolution and subsequent ion formation in a broad range of complex matrices at low sample concentrations. Additionally, these findings contribute to the progress of biomodelling and dosimetry approaches to close the gap between in vivo and in vitro results. Thus, the results deepen the understanding of NP toxicity.

Chemicals and Materials
If not noted, otherwise stated chemicals were purchased at Merck (Darmstadt, Germany) with at least 98% purity. The sample handling was performed in 50 mL tubes (Corning 430828, Corning, New York, NY, USA). Stock solutions of copper oxide NPs (544868, Merck, Darmstadt, Germany) were prepared in accordance with the Nanogenotox protocol: "Final protocol for producing suitable manufactured NMs exposure media" (October 2011). Shortly, a 2.56 mg mL −1 stock dispersion was freshly prepared by pre-wetting the powder with 0.5% (v/v) ethanol (96%) followed by the addition of Milli-Q water (Milli-Pore gradient, Merck, Darmstadt, Germany) containing 0.05% BSA. The dispersion was sonicated for 5 min and 9 s at an amplitude of 10% with a probe sonifier (200 W Bandelin Sonopuls HD 2200, BANDELIN Electronic GmbH & Co. KG, Berlin, Germany). The sample was cooled in an ice-water bath during sonication [44].
For the enhanced Gamble's solution, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (16:0 PC; DPPC) was provided by Avanti Polar Lipids, Inc. (850355C via Merck, Darmstadt, Germany). To ensure a fine dispersion of DPPC, a sonication was performed as described in the Nanongenotox protocol after its addition. The different media were prepared according to the literature [41]. As a commonly used cell culture medium, RPMI 1640 (Merck, Darmstadt, Germany) was used. Similar to a normal use case, 10% fetal bovine serum and 1% antibiotics (10,000 µg mL −1 streptomycin and 10,000 units mL −1 penicillin) were added to the RPMI. The freshly prepared NP stock solutions were further diluted with regard to the different techniques. Final concentrations of 25 mg L −1 , 20 µg L −1 and 1 µg L −1 in the media were used for the DLS measurements, for the microwave digestions and for the direct measurements, respectively. Standard solutions of dissolved Au and Cu (1000 mg L −1 , TraceCERT, Sigma-Aldrich, Darmstadt, Germany), 3.5% nitric acid (technical grade 70% v/v, VWR, Darmstadt, Germany) were purified in a douPur quartz sub-boiling point apparatus (MLS GmbH, Leutkirch im Allgäu, Germany), and the other media were used for the calibration solutions (1, 2, 5 and 10 µg L −1 ). The reference Au-NPs NIST 8012 (30 nm, Gaithersburg, MD, USA) were diluted in the respective media to concentrations down to 50 ng L −1 to determine matrix-matched transport efficiencies of the DI spICP-MS set-up.

Experimental Details
All samples were stored at 20°C. At defined time points (approx. at 0, 1, 2, 4, 6, 8, 24, 48 and 168 h) samples of CuO NP dispersions inside the different media were characterized with DLS after vortexing the samples.
The sample preparation for the acid-assisted microwave digestions consisted of a vortexing step, after which a volume of 1 mL was extracted into an Eppendorf tube.
Afterwards, the samples were centrifuged at 200 rpm for 5 min to force the sedimentation of the NPs. Directly following this step, the top 500 µL of the samples were separated, and the pellet inside the leftover solution was redispersed through vortexing. As a next step, the solutions were digested by adding 69% HNO 3 and 30% H 2 O 2 and processing them at 200°C and 160 bars. Afterward, a 1:50 dilution with Milli-Q water was performed prior to the ICP-MS measurements.
In the case of DI spICP-MS, measurements were performed more frequently (between 15 min to 1 h) to a total of 20-30 time points per experiment. Each experiment was performed with three independent replicates and additional media blanks.

Instrumentation and Data Evaluation
The DLS analysis was performed with a Zetasizer from Malvern Panalytics (Kassel, Germany). Thermal equilibration time was set to 60 s at 25°C. Each value and time point represents the average of five individual DLS measurements using automatic optimization of analytical conditions and data treatment by general-purpose size analysis. NP sizes were determined as z-average and intensity mean. Furthermore, the polydispersity index (PDI) and the count rate were visualized. Bar plots of these values average the measurements below 12 h against the determined averages after several days (i.e., 1-7). Therefore, the relative standard deviation (RSD) is influenced by the uncertainty of the individual measurements and the processes over time.
For the evaluation of ionic contents of the digested solutions, measurements were performed with a quadrupole ICP mass spectrometer (iCAP Q, Thermo Fisher Scientific GmbH, Dreieich, Germany) equipped with a PFA ST Nebulizer, a quartz cyclonic spray chamber and a 2.5 mm quartz injector (all from Thermo Fisher Scientific). The gas flow for the cool gas (Ar) and the auxiliary gas (Ar) were set to 14 and 0.65 L min −1 , respectively. The sample flow rate was determined to be 0.34 mL min −1 . 59 Co and 65 Cu were analyzed using the collision cell technique at 5 mL min −1 (collision gas flow with 93% He and 7% H 2 ).
Measurements of the CuO NPs were performed with the 65 Cu signal to avoid possible interferences from salts. Due to the complete dissolution of NPs as part of the microwave digestion, the ionic contents of the supernatant and the sedimented NPs were evaluated. These values were averaged to the total ionic content, and the particle-based part was calculated by taking into account the measured ionic concentration of the supernatant. In addition to the bar plots with long and short-term values, the ionic distributions at the individual time points were visualized to highlight the different influences of the media on the NP dissolution.
For single particle analysis of the NP solutions, a quadrupole ICP-MS (Thermo Scientific XSERIES II, Thermo Fisher Scientific, Waltham, MA, USA) with a PFA ST Nebulizer, a quartz cyclonic spray chamber and a 2.5 mm quartz injector (all from ESI Elemental Service and Instruments GmbH, Mainz, Germany) were used. Using the time-resolved analysis mode for data acquisition, intensities were collected as a function of time (counts per dwell-time interval). The acquisition time for each run was set to 70 s with a dwell time (or data acquisition rate) of 3 ms. The gas flows for the plasma, the nebulizer and the auxiliary (all Ar) were set to 13, 0.89 and 0.7 L min −1 , respectively. All measurements were performed using the collision cell technique to avoid polyatomic interferences with the same collision gas and flow rates as mentioned before. The injection system was a modified microFAST MC system (ESI Elemental Service and Instruments GmbH, Mainz, Germany) with a single loop (300 µL) set-up and 2 independent syringes (5 mL and 500 µL). A sample loop of 300 µL was used with a flow rate of 500 µL min −1 to ensure a transient signal independent from the media composition. The measurement scripts were adjusted in accordance with the set-up characteristics and to reduce the overall measurement time. Routine tuning of this system was performed with an adapted routine and the commonly used Tune F solution (ESI Elemental Service and Instruments GmbH, Mainz, Germany) to maximize the signal intensity of 59 Co. Results of an aqueous 5 µg L −1 Au solution, as described above, were used to evaluate the performance at different flow rates (50, 75, 100, 150, 300 and 500 µL min −1 ) based on the averaged duration, intensity and its RSD.
Further experiments were performed with a flow rate of 500 µL min −1 to avoid the clogging effects of higher viscosity media. For the evaluation of these parameters and the determination of matrix-matched transport efficiencies, the 197 Au signal was observed. The determination of the threshold between ionic background and particle signals and the respective counting of the observed particle number were performed systematically, based on the corrected average intensities of the ionic background in an automated fashion. First, the transient signal of the individual data set was divided into 10 consecutive pieces (typically 1000 data points each) to take drifts during a single measurement into account. Next, the total average value and average of each piece was determined.
To determine the threshold, either 3 or 5 times the RSD (σ) was added to each average [45,46]. The individual point values above their respective limit were taken out and the whole procedure was repeated 5 times. The final limit of each timepiece was used to determine the particle number of each data piece, respectively. The averaged particle number per measurement was used to determine the matrix-matched transport efficiency.
Similar to the automated approach described above, the particle numbers, as well as the particle size and ionic background, were calculated at each time point. The basis of this approach is the excel spreadsheet from RIKILT Wageningen UR available at https: //www.wur.nl/en/show/Single-Particle-Calculation-tool.htm (accessed on 21 September 2021). To obtain the total copper content for the visualization similar to the microwave digestions, the ionic content and particulate content were summed up. The particle numbers were corrected using the matrix-specific transport efficiencies, and the concentration was derived based on the actual volume used.

Nanomaterial Characterization of Hydrodynamic Parameters in Complex Media
An overview of different observed key parameters of DLS measurements of complex media spiked with dispersed CuO NP determined at different time points is presented in Figure 1.

Results of dynamic light scattering (DLS) measurements of CuO nanoparticle (NP)
dispersions in several media at different time points  The observed size reflects only the hydrodynamic diameter of the particles inside the samples. Therefore, the observed size corresponds to the overall size of a particle, e.g., their metal/metal oxide and their so-called "corona" formed by matrix components. Furthermore, scattering particles influence the overall results with different contributions based on their size. This is, to a certain extent, reflected in the comparison of the sizes obtained by the z-average compared to the intensity mean-based hydrodynamic radii (see Figure 1). The polydispersity index (PDI) provides a general idea of the complexity of the size distribution of the sample. A PDI value below 0.05 represents a monodisperse sample, while a PDI above 0.7 reflects a broad size distribution of particle characteristics of a polydisperse sample. In combination with the count rate of the respective averaged measurements, these two parameters can be used to evaluate the measurement conditions in general and interpret the obtained results. Changes in the count rate provide a hint about the stability of the particles inside the solution, despite its automatic optimization. An average of the results based on the time points at 0, 1, 2, 4, 6, 8 h (see Figure 1, orange bars) and the results obtained after 1, 2 and 7 days (see Figure 1, green bars) highlights the different behavior of the same amount of CuO NPs inside the different media. The principle behind the size determination leads to challenges of measurements with ALF, GS and EGS matrices. The influence of newly formed or modified particulate content can be determined by significantly higher PDI values and more fluctuating count rates between individual repeats and time points. Another point, which has to be taken into account, is the less stabilizing effect of these matrices, which contributes to the higher standard deviations. These findings indicate changes in the particle content of the GS, ALF and EGS solutions. Due to insufficient particle numbers, the measurements of the blank media could not be used as a reference for background levels of particulate matter formed out of the media itself. In the case of GS, it was reported that the ionic content is not stable over time. This might lead to particle formation and agglomeration of the ionic content. In case of the EGS solution, the decreasing PDI and counts, in combination with the increasing size (intensity-based), can be attributed to the formation of lipid structures, such as liposomes, inside the solution or the formation of bigger NP-lipid clusters. A similar effect could be observed in experiments with gold NPs and liposomes [43]. These agglomerates may be formed over time and enhance the NP agglomeration too. Lastly, the results of the ALF solutions can be explained with a loss of particles derived from the significant decrease in counts between the time points. The pH-dependent dissolution of the CuO NPs in this acidic media points in this direction. This can be observed in the determined sizes based on the intensity mean distributions too. In comparison, the solutions containing water, BSA or RPMI cell culture medium provide similar responses for the size determination approaches with slight changes over time. The PDI shows slight differences between these media, as assumed from their composition. The count rates provide a stable measurement with smaller changes over time. A conclusion might be that the stabilizing effect of the bio-molecules and the buffer solution lead to better working conditions in terms of DLS measurements. With this summary of the results, the advantages of DLS are highlighted, which are to provide a first and fast insight into changes or their absence. A remaining challenge is correlating the observed parameters correctly to the underlying mechanisms. This is especially true for complex samples and matrices, as shown here.

Microwave Digestions of CuO NP Dispersions in Different Media
A different approach to dealing with complex media and their potential limiting effects is the separation of NP sample material from dissolved ionic content by centrifugation, the reduction in the matrix by microwave digestion and a successive quantification by ICP-MS of the respective ionic contents. This reduces the burden on the characterization techniques at the cost of a loss of information. To avoid uncertainties due to the additional preparation steps, the individual absolute ionic contents were normalized to the respective total content values. The big advantage of this strategy is its high versatility and robustness. The results presented in Figure 2 show an increase in ionic content in all media, albeit only slightly in water.
Ionic content of CuO NP dispersions in several media after centrifugation and microwave digestion at different time points Accordingly, the particulate content decreases after 8 h significantly. Due to the sample preparation approach in combination with the ICP-MS characterization, these results reflect different processes compared to the DLS measurement. The internal comparison and interpretation of the time-resolved data allow for a more detailed look into the individual measurements of the different media (see Figure 3). This is possible due to the lower individual uncertainties at the respective time points compared to the DLS measurements.
Ionic copper content of CuO NP dispersion at multiple time points in different media after their centrifugation and microwave digestions The trends observed with the DLS measurements are only partly visible due to the different parameters observed. The water solution shows high uncertainties over all time points, which prevents an in-depth interpretation. One point, which needs to be considered, is the additional uncertainty originating from the increased number of working steps. This may lead to a bigger number of experimental errors. Furthermore, the centrifugation-based separation of dissolved and particulate content depends on the density of the NPs and takes time. Both considerations might slightly influence the result toward higher ionic values. However, a slight loss of particulate content can be seen relative to the ionic background values. In contrast, the RPMI-containing samples maintain a stable ratio for several hours after the particles are dispersed homogeneously. A similar pattern can be more clearly observed in the EGS medium. A possible explanation is the stabilizing effect of bigger molecules from the media itself, e.g., proteins and lipids. This is in line with the exclusion of pH effects due to the buffered solutions. In terms of the GS solution, which only contains smaller ionic molecules, a steady increase in ionic content can be observed. Lastly, the fastest dissolution can be seen in ALF. This further verifies the pH-driven dissolution of the CuO NPs in this acidic medium. Due to the reduced matrix influences caused by the microwave digestion, this approach provides a robust way to answer questions about the dissolution of metal NPs in complex media. However, if the characterization of the individual NPs (e.g., size (distribution), shape, number) is required, this method cannot supply sufficient information.

Development of a Direct-Injection Single Particle Inductively-Coupled Plasma Mass Spectrometry System
A method that combines the high sensitivity of the ICP-MS with information about the individual NPs is the spICP-MS. This approach is typically performed in diluted aqueous samples with diluted nitric acid as a carrier to enhance the overall performance. The challenge of this method is based on carry-over effects of the diluted acid due to the autosampler system, as well as the higher necessary sample concentration. These challenges were dealt with with the DI spICP-MS approach and pure water as the carrier solution. First, the feasibility of the DI spICP-MS set-up, despite the different complex matrices and the loss in overall sensitivity due to a non-acidic measurement environment, is evaluated and addressed. This includes a general optimization of the overall set-up and experimental parameters with ionic and particulate standards. Following these experiments, CuO NP dispersions in the complex media were characterized.
For these measurements, a modified spICP-MS set-up with 2 syringes (5 mL and 500 µL) for exact volume handling, stable performance and fast measurements were used. The challenge of such a direct injection system is the effect of the different undiluted media on the experimental set-up and process parameters. Additionally, matrix effects during the transport, spray, nebulization and ionization processes might further influence the results. Therefore, the goal was a reduction in the necessary sample volume to reduce the burden on the instrument. After the general set-up was configured, further enhancement of signal intensities was achieved due to small amounts of air at the start and end of the sample volume inside the sample loop. With these arrangements, transient signal profiles could be achieved with flow rates between 50 and 500 µL per minute (see Figure 4). The general effects observed with an increasing flow rate are a smaller time frame of the transient signal with an increased signal intensity and a decreased relative standard deviation. In this example, a plateau of those effects is reached above a flow rate of roughly 200 µL/min. However, the viscosity of the different media might inhibit the application of lower flow rates due to clogging of the nebulizer. To have a fast, robust and comparable set of parameters, a flow rate of 500 µL/min was chosen. The downside of this approach is the relative loss of signal intensity compared to the average measurement duration. In theory, this limits the technique to lower NP number concentrations to avoid possible double NP events. A dilution of the sample would solve these issues and reduce the analytical challenges.

Exemplary measurements of ionic gold solution at different flow rates
and key parameters of these experiments A general overview of the results of NP measurements (Au NIST 8012) in the individual media is presented in Table 1. Another significant challenge of spICP-MS is the selection of a limit for the ionic background. This is especially true in the case of a naturally occurring ionic content, dissolving particulate content or small particle sizes, which leads to low spike intensities. To avoid the influence of a subjective determination of this limit, an automated approach was performed. The limit is based on the corrected ionic average and a safety margin of 3 σ or 5 σ RSDs of the individual time frame to calculate the number of NP events. As shown in Table 1, successful measurements and data evaluation could be performed in all media. Due to this mathematical approach, a systematic error can be observed when sizes are compared to the certified size of the particles (approx. 30 nm). This is especially true for the higher safety margin of the 5 σ approach. Therefore, the 3 σ calculations were chosen for the data interpretation of the CuO NP experiments.
Another potential systematic error might be due to ionic background values below the electrical noise. In these cases, an overestimation of the number of events due to these noise values is to be expected. Nonetheless, the observed particle numbers per measurement and general considerations are in line with the recommendation of the ISO/TS 19590:2017 [43].
In conclusion, for a quantitative evaluation, careful calibration and reference measurements are necessary to compensate for these effects. Another option is to further optimize the experimental set-up and parameters for each media, respectively. The 3 σ or 5 σ approach to determine the ionic background levels individually for each time frame and measurement grants an objective way to handle complex cases in an automated, reproducible way. This advantage becomes more pronounced when comparing results over time. Then, the day-today performance of the instrument might influence the measurement results significantly, in addition to the present matrix effects.

Nanoparticle Characterization in Complex Media with DI spICP-MS
To evaluate the presented approach (see Section 3.3) of an objectively evaluated data treatment procedure in combination with a direct-injection measurement method, experiments were performed with the CuO NPs dispersed in a broader variety of media. A summary of these experiments is shown in Table 2 and Figures 5 and 6. In all different media, successful measurements could be carried out. Exemplary time-resolved results of the respective media are provided in the supplementary materials (see Figures S1-S6). Due to the reduced time for each individual measurement, the lower required volumes and the automatized procedure, a better time resolution could be achieved. However, the matrix influences the sensitivity of the measurements, the observed ionic background levels, as well as the calculated transport efficiencies. A comparison with bar plots shows the advantages of this technique (see Figure 5). Results of direct-injection single particle (DI sp)ICP-MS measurements of CuO NP dispersions of several media at different time points The size of the particles based on the time points at the start and after 7 days provides a significant difference in BSA, buffered systems, ALF and GS media. In the case of BSA and RPMI, the number of particles determined after 7 days was below the recommended limit of the ISO/TS 19590:2017; therefore, these values are written in parentheses in Figure 5 [43]. It is important to note that in the ICP-MS-based evaluation, the element-specific information is the basis of the size information. This explains different findings compared to the DLSdetermined hydrodynamic sizes in the artificial lung surfactant media. For the pure water, EGS, and cell culture media, no significant change in size could be determined. This is in line with the findings of the DLS measurements. Although, the sizes determined with the DI spICP-MS approach are influenced by matrix effects (see Figure 5, RPMI), which might be attributed to the copper content of the media. The blank media measurements revealed the RPMI as the origin of dissolved copper, leading to a high ionic background. The observed sizes are in most cases significantly higher than the size provided by the manufacturer (<50 nm). Possible reasons are related to the influence of the different media and the size distribution of the CuO NPs. On the one hand, the particle agglomeration itself can be influenced by the constituents of the media. Additionally, matrix effects lead to a reduced sensitivity of the DI spICP-MS. Therefore, smaller-sized NPs might not be detected, and thus, the calculated particle size based on the detected bigger NPs is above the value claimed by the manufacturer. Furthermore, similar to the DLS measurements, this approach provides more information than a calculated spherical size (see Table 2 and Figure 5b,c). The comparison of the ionic content at the two time points shows a stable amount in several dispersions over the week. Nonetheless, the matrix effects strongly influence the sensitivity, which leads to high uncertainties for the DI spICP-MS. In contrast, ALF shows a clear dissolution, with a significantly increased ionic copper content. This behavior is only partly reflected in the observed count rates and PDI of the DLS data (see Figure 1c,d). Finally, the visualization of the particle number concentrations highlights an NP decrease over time.
In combination with the determined ionic content of the individual time points, losses can be attributed to dissolution or other unidentified processes In the case of EGS and RPMI, the matrix influences reduce this effect. A more detailed look into the dissolution process was provided with the microwave digestion approach (see Figure 3).
A comparison of these findings with the DI spICP-MS results is provided in Figure 6. Due to the calculation of ionic content based on the determined particle numbers, a significant increase in ionic content attributed to ionic background values can be observed in all media used in the microwave-assisted approach. However, these findings are influenced by matrix effects such as higher background values and lower sensitivity, which cause a higher standard deviation. These effects are reduced due to the digestion and dissolution in the conventional ICP-MS experiments. Nonetheless, the additional information provided by the DI spICP-MS approach, in combination with the reduced efforts, lead to unique applications and insights due to characterized particles inside the media. This can be seen in the right-hand-side graphs of Figure 6, where a visualization of the determined averaged particle numbers at several time points shows different types of kinetics inside the various media. To sum up and generalize these results: the DI spICP-MS approach proves to be a fast and versatile technique to observe metallic NPs directly. This technique, with the systematic data evaluation strategy in combination with the fast measurements, offers a strong tool to determine particle fate. Lastly, detailed optimizations of the individual process parameters can increase the performance further for a certain analytical question (e.g., a specific medium, ionic concentration, particle number or distribution) as shown in Figure 4.

Conclusions
A key point of this study is the discussion of different techniques and their applicability for different analytical questions with relation to NP characterization. Each technique has its individual advantages and challenges. It is, therefore, important to define which parameters must be investigated at which measurement conditions, including the amount of time and resources available. The analytical question itself includes the differentiation between quantitative and qualitative information and, additionally, the requirements in terms of precision and reproducibility. As a baseline, pure aqueous solution medium was selected. To increase the analytical challenges step-wise, BSA-containing solution and sodium hydrocarbonate buffer were chosen. Due to their application in toxicological experiments, a representative cell culture medium was selected. With regard to experimental dissolution scenarios and their biological relevance, ALF, GS and EGS were used. The observed time frame of 7 days was sufficient to characterize NP behavior. The results of the CuO NP dispersions in pure aqueous media showed a significant loss of particles. The determined sizes and ionic contents over time were in the same range due to the high uncertainty of the measurements. Therefore, the origin of this particle loss cannot be determined. The experiments performed with NP dispersions in BSA-containing media demonstrate a clear loss of NP numbers as well. Additionally, the ionic content increased significantly over time, and the determined NP size decreased. Therefore, the assumption can be made that NP dissolution and agglomeration of bigger particles followed by with their sedimentation took place. Other processes leading to particle loss might have occurred in all media too. The NP dispersions with buffered media showed a less pronounced loss of particles. The observed NP sizes and the ionic content increased slightly. This indicated a slower dissolution due to the avoidance of a pH-driven dissolution mechanism. Furthermore, a significant loss of bigger particles, decreasing the overall particle size, can be excluded. The CuO NP dispersions with the acidic ALF showed a fast loss of particles, a strong increase in ionic content and a slight decrease in particle size. A fitting explanation was the existence of a pH-driven dissolution. However, agglomeration and sedimentation were likely to happen as well. The experiments with the GS media indicated a significant loss of particles, a decrease in size and an increase in ionic content, as shown by the MW ICP-MS data. Similar to the conclusions of the BSA-containing NP dispersions, dissolution and a combination of agglomeration and sedimentation were likely to take place. The experiments of NP dispersion in EGS showed only a slight loss of particles and a decrease in size, while the MW ICP-MS results revealed an increase in ionic content. Here, dissolution took place, but the agglomeration behavior is mitigated compared to the observed behavior of the experiments with the GS media. Lastly, the CuO NP dispersions in RPMI media showed the lowest starting number of particles, and due to the uncertainties, no loss can be determined. However, the digestion experiments revealed a slight dissolution. Due to the same observed sizes, sedimentation processes seem unlikely. An overview of these mentioned aspects and their impact on the applicability of the shown techniques is provided in Table 3. Additional NP traits, which might impact the analysis, such as shape, were neglected here as no further data regarding these topics were collected. Based on the results presented here, it can be concluded that the DLS measurements act as a fast, low-cost method, albeit unknown matrix influences might be challenging. The microwave digestion proves to be a robust and reliable method at the cost of information loss and time-consuming sample preparation steps. The newly introduced combination of a DI spICP-MS approach with an automatized data evaluation procedure produces fast and reproducible results for a broad range of complex media. Nevertheless, matrices can influence the overall performance and instrumental strain. Therefore, an in-depth optimization must be performed for a quantitative characterization. Further technical improvements or adapted configurations may lead to a higher precision or multi-elemental analysis. This approach is especially promising in differentiating between ionic-based observed toxicity effects of metallic NPs and particulate-driven processes. This is exemplarily shown in the acidic scenario of the ALF compared to the other matrices. This new information could expand the knowledge about ongoing NP-related processes in different environments and allows the development and utilization of less toxic NP.

Data Availability Statement:
The data presented in this study are available upon justified request.